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ABSTRACT 

Three different methods of measuring cosmology with gamma-ray bursts 
(GRBs) have been proposed since a relation between the 7-ray energy E 1 of 
a GRB jet and the peak energy E p of the vF v spectrum in the burst frame was 
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already best fitted for this cosmology is considered. We apply this method to a 



sample of 17 GRBs, and obtain the mass density Q M = 0.15±g;^ (1 a) for a flat 
ACDM universe. In Method II, to calculate the probability for some certain cos- 
mology, contributions of all the possible E 1 — E p relations that are best fitted for 
their corresponding cosmologies are taken into account. With this method, we 
find a constraint on the mass density 0.14 < Qm < 0.69 (la) for a flat universe. 
In Method III, to obtain the probability for some cosmology, contributions of all 
the possible E y — E p relations associated with their unequal weights are consid- 
ered. With this method, we obtain an inspiring constraint on the mass density 
0.16 < tt M < 0.45 (la) for a flat universe, and a xLf = 19.08/15 = 1.27 for the 
concordance model of Qm = 0.27. Compared with the previous two methods, 
Method III makes the observed 17 GRBs place much more stringent confidence 
intervals at the same confidence levels. Furthermore, we perform a Monte Carlo 
simulation and use a larger sample to investigate the cosmographic capabilities 
of GRBs with different methods. We find that, a larger GRB sample could be 
used to effectively measure cosmology, no matter whether the E 1 — E p relation 
is calibrated by low- z bursts or not. Ongoing observations on GRBs in the Swift 
era are expected to make the cosmological utility of GRBs progress from its 
babyhood into childhood. 



Subject headings: gamma rays: bursts — cosmology: observations — cosmology: 
distance scale 
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1. Introduction 

The traditional cosmology has been revolutionized by modern sophisticated observation 
techniques in distant Type la supernovae (SNe la) (e.g. Riess et al. 1998; Schmidt et 
al. 1998; Perlmutter et al. 1999), cosmic microwave background (CMB) fluctuations (e.g. 
Bennett et al. 2003; Spergel et al. 2003), and large-scale structure (LSS) (e.g. Allen et al. 
2003; Tegmark et al. 2004). Each type of cosmological data trends to play an unique role 
in measuring cosmology. In modern cosmology, it has been convincingly suggested that the 
global mass-energy budget of the universe, and thus its dynamics, is dominated by a dark 
energy component, and that the currently accelerating universe has once been decelerating 
(e.g., Riess et al. 2004). The cosmography and the nature of dark energy as well as its 
evolution with redshift are one of the most important issues in physics and astronomy today. 

Gamma-ray bursts (GRBs) are the most intense explosions observed so far. They are 
believed to be detectable up to a very high redshift (Lamb & Reichart 2000; Ciardi & Loeb 
2000; Bromm & Loeb 2002; Gou et al. 2004), and their high energy photons are almost 
immune to dust extinction. These advantages would make GRBs an attractive cosmic probe. 

From the isotropic-equivalent peak luminosity L iso — variability (or spectral lag) relation 
(Fenimore & Ramirez- Ruiz 2000; Norris, Marani, & Bonnell 2000), the standard energy 
reservior E 1 of GRB jets (Frail et al. 2001), the L iso — peak energy E p of the vF v spectrum 
in the burst frame relation (Lloyd-Ronning & Petrosian 2002; Lloyd-Ronning & Ramirez- 
Ruiz 2002; Yonetoku et al. 2004), the isotropic-equivalent energy E iso — E p relation (Amati 
et al. 2002) to the beaming-corrected energy E y — E p relation (Ghirlanda et al. 2004a, 
Ghirlanda relation hereafter), GRBs are towards more and more standardizable candles. 
However, these relations in GRBs have not been calibrated by a low- z GRB sample, so one 
should look for a method which is different from the "Classical Hubble Diagram" method in 
SNe la. 

The luminosity relations with the variability and spectral lag make GRBs a distance 
indicator in the same sense as Cepheids and SNe la, in which an observed light-curve property 
can yield an apparent distance modulus (DM). Schaefer (2003) considered these two relations 
for nine bursts with known redshifts and advocated a new cosmographic method (hereafter 
Method I) for GRBs. In Method I, one first calibrates the two relations with the observed 
sample for a certain cosmology, and then applies the best-fit relations back to the observed 
sample to obtain a x 2 or a probability P oc exp(— % 2 /2) for this cosmology. Similar to the 
brightness of SNe la, the energy reserviors in GRB jets are also clustered, but they are not 
fine enough for precise cosmology (Bloom et al. 2003). Amati et al. (2002) found E iso oc E p k 
(k ~ 2) relation from 12 BeppoSAX bursts. The HETE-2 observations confirm this relation 
and extend it to X-ray flashes (Sakamoto et al. 2004a; Lamb et al. 2004). In addition, it 
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also holds within a GRB (Liang, Dai & Wu 2004). The Ghirlanda relation is written as 
(£' 7 /10 50 ergs) = C(-E p /100keV) a , where a and C are dimensionless parameters. Theoretical 
explanations of this relation include the standard synchrotron mechanism in relativistic 
shocks (Zhang & Meszaros 2002; Dai & Lu 2002) together with the afterglow jet model, 
or the emission from off-axis relativistic jets (Yamazaki, Ioka & Nakamura 2004; Eichler & 
Levinson 2004; Levinson & Eichler 2005). This relation could also be understood due to 
comptonization of the thermal radiation flux that is advected from the base of an outflow in 
the dissipative photosphere model (e.g. Rees & Meszaros 2005). If these explanations are 
true, the Ghirlanda relation appears to be intrinsic. Thus, Dai, Liang & Xu (2004; DLX04) 
considered the Ghirlanda relation for 12 bursts and proposed another cosmographic method 
(hereafter Method II) for GRBs. In Method II, one makes marginalizations over the unknown 
parameters in the Ghirlanda relation to obtain a f or a probability P oc exp(— x 2 /2) for 
a certain cosmology. Following Schaefer's method, Ghirlanda et al. (2004b; GGLF04) and 
Friedman & Bloom (2004; FB04) also investigated the same issue but used different GRB 
data. Recently, Firmani et al. (2005; FGGA05) considered the Ghirlanda relation for 15 
bursts and proposed a Bayesian approach for the cosmological use (hereafter Method III). 
In Method III, to obtain the probability for a certain cosmology, one considers contributions 
of all the possible E 1 — E p relations associated with their unequal weights. The detailed 
procedures of the three methods are shown in §2.1, which indicate that Method III is the 
optimized one. 

As analyzed previously, due to the lack of low- z GRBs, Methods I, II and III are different 
from the "Classic Hubble Diagram" method in SNe la. In this paper, we investigate the 
constraints on cosmological parameters from the observed 17 GRBs with different methods. 
Because the present GRB sample is a small one, it is necessary to use a large simulated 
sample, which may be established in the Swift era, to discuss the cosmographic capabilities 
with different methods. 

This paper is arranged as follows. In §2, we describe our analytical methods and data. 
The results from the observed GRB sample are presented in §3. In §4, we perform Monte 
Carlo simulations and analyze the results from the simulated GRB sample. Conclusions and 
discussion are presented in §5. 
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2. Method and Sample Analysis 

2.1. Method Analysis 

According to the relativistic fireball model, the emission from a spherically expanding 
shell and from a jet is similar to each other, if the observer is along the jet's axis and the 
Lorentz factor of the fireball is larger than the inverse of the jet's half-opening angle 9; 
but when the Lorentz factor drops below the jet's afterglow light curve is expected to 
present a break because of the edge effect and the laterally spreading effect (Rhoads 1999; 
Sari, Piran & Halpern 1999). Therefore, together with the assumptions of the initial fireball 
emitting a constant fraction r^ 7 of its kinetic energy into the prompt 7-rays and a constant 
circumburst particle density n, the jet's half-opening angle is derived to be 
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where -E iso ,52 = -E iS o/10 52 ergs, t ijd = ij/lday, n = n/lcm" 
equivalent 7-ray energy of a GRB is given by 
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where S* 7 is the fluence (in units of ergcm~ 2 ) received in an observed bandpass and the 
quantity A; is a multiplicative correction of order unity relating the observed bandpass to a 
standard rest-frame bandpass (1-10 4 keV in this paper) (Bloom, Frail & Sari 2001). The 
energy release of a GRB jet is thus given by 



E 1 = (1 - cos#)i?i S o, 
with the fractional uncertainty (FB04) being 
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The Ghirlanda relation is 
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where a and C are assumed to have no covariance, and E p = E° hs (l + z). Combining 
Eqs. (1), (2), (3) and (6), we derive the apparent luminosity distance with the small angle 
approximation (i.e., 6 <C l) 1 as 
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Assuming that all the observables are independent of each other and their errors satisfy Gaus- 
sian distributions, we derive the fractional uncertainty of the apparent luminosity distance 
without the small angle approximation, 
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For simplicity, we consider i] y = 0.2 and a Vl = throughout this paper (Frail et al. 2001). 
The apparent DM of a burst can be given by 



with the uncertainty of 



/i obs = 5\ogd L + 25, 
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On the other hand, the theoretical luminosity distance in A-models (Carroll, Press & 
Turner 1992) is given by 
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where = 1 — — &a, and "sinn" is sinh for > and sin for £l k < 0. For £l k = 0, 
equation (11) degenerates to be c(l + z)!!^ 1 times the integral. 



1 This approximation is valid because |(1 — cosd — 9 2 /2)/(l — cos#)| < 1% when 9 < 0.35 rad, and < 0.4% 
when 9 < 0.22 rad. 
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Usually, the likelihood for the parameters VLm and Q\ can be determined from a x 2 
statistic, where 



where the dimensionless Hubble constant h = ifo/100 kms _1 Mpc _1 is taken as 0.71. If the 
Ghirlanda relation could be calibrated by low- z bursts, the above x 2 statistic becomes the 
same as that in SNe la, that is, 



where h should be marginalized. 

The procedures of Method I , II and III are as follows: 
Method I (see S03, GGLF04, FB04) 

The procedure of this method is to (1) fix fij = (^m,^a)i, (2) calculate fi^ and E 1 
for each burst for that cosmology, (3) best fit the E 1 — E p relation to yield (a, C)j and 
(a a /a, ac/C)i, (4) substitute (a, C)i and (o" a /a, cr c /C)i into Eqs. (7) and (8), and thus de- 
rive /i bs and o" Atobs for each burst for cosmology Qi, (5) calculate x 2 for cosmology Qi by com- 
paring /i th with /x bs, °"/i obs5 and then convert it to the probability by P(fij) oc exp(— x 2 (^«)/2) 
(Riess et al. 1998), (6) repeat Steps 1—5 from % — 1 to % — N to obtain the probability for 
each cosmology. Therefore, Method I is formulizd by 



Method II (see DLX04) 

The procedure of this method is to (1) fix (2) calculate fi th and for each burst 
for that cosmology, (3) best fit the _E 7 — E p relation to yield (a,C)j and (a a /a, a c /C)i, (4) 
substitute (a, C)i and (cr a /a, ac/C)i into Eqs. (7) and (8), and thus derive /i Q b s and <r Mobs 
for each burst for cosmology Qi, (5) repeat Steps 1—4 from i — 1 to i — N to obtain /i t h, 
/i bs and <7 Mobs for each burst for each cosmology; (6) re-fix Qj, (7) calculate x 2 (r2 J |f2 i ) by 
comparing fi t ^(Qj) with /i b s (^j), a^^Qi), and then convert it to a conditional probability 
by P(Qj\Qi) oc exp(— x 2 (fij |r2j)/2), (8) repeat Step 7 from % — 1 to % — N to obtain the 
probability for cosmology Qj by P{Qj) oc J^exp(— x 2 (fij|fij)/2), (9) repeat Steps 6—8 from 



j — 1 to j — N to obtain the probability for each cosmology. Method II is described by 
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P{Qi) = P(Qil^) (i — 1, N). 
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Method III (see FGGA05) 

Method III is an improvement of Method II. Its key idea is to consider unequal weights 
for different E 1 — E p relations, i.e., unequal weights for different conditional probabilities 
P(Qj\Qi). Therefore, the first seven steps of this method are the same as those of Method 
II. The follow-up procedure is to (8) repeat Step 7 from i — 1 to i — N to obtain an iterative 
probability for cosmology Qj by P itc (f^) oc ex P( _ X 2 (m^i)/2) x P ini (fij) (here the initial 

i 

probability P im (f2) for each cosmology is regarded as equal; e.g. P im (f2) = 1), (9) repeat 
Steps 6—8 from j — 1 to j — N to obtain an iterative probability P lte (f2) for each cosmology; 
(10) replace P im (f2) on Step 8 with P ltc (f2) on Step 9, then repeat Steps 8—9, and thus reach 
another set of iterative probabilities for each cosmology, (11) run the above cycle again and 
again until probability for each cosmology converges, i.e. P lte (f2) =>- P fin (f2) after tens of 
cycles. 

In this method, to calculate the probability P fin (f2j) for a favored cosmology, we con- 
sider contributions of all the possible E 1 — E p relations associated with their weights. The 
conditional probability P(flj\Qi) denotes the contribution of some certain relation, and 
P (fli) weights the likelihood of this relation for its corresponding cosmology. Therefore, 
the Bayesian approach can be formulized by 

N 

p Ba {n j ) = Y,P(to j \n i )xP Ga (n i ) (j = i,n). (16) 

i=i 

However, FGGA05 took different calculations for the conditional probability P(Qj\Qi). 
Making use of the incomplete gamma function, they transformed x 2 {^j\^i) m t° hs cor- 
responding conditional probability P(Qj\Qi). Actually, once the parameters (a,C)i and 
(cr a /a, a c /C)i of the Ghirlanda relation are calibrated for cosmology Jlj, they become 
"known" for the cosmic model Qj. So herein the meaning of x 2 (^j\^i) is the same as that 
in SNe la. In this paper, we redefine the conditional probability P(Qj\Qi) by the formula of 
oc exp(-x 2 (^-|^)/2). 

2.2. Sample Analysis 

The great diversity in GRB phenomena suggests that the GRB population may consist 
of substantially different subclasses (e.g. MacFadyen & Woosley 1999; Bloom et al. 2003; 
Sazonov et al. 2004; Soderberg et al. 2004). To make GRBs a standard candle, a homogenous 
GRB sample is required. The most prominent observational evidence for a GRB jet is its 
temporal break in their afterglow light curves. For some bursts, e.g., GRB030329, their 
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temporal breaks are observed in both optical and radio bands. Berger et al. (2003) argued 
that these two breaks are caused by the narrow component and wide component of the jet in 
this burst, respectively, indicating that the physical origins of the breaks in the optical band 
and in the radio band are different. In addition, the radio afterglow light curves fluctuate 
significantly. For example, in the case of GRB970508, the light curve of its radio afterglow 
does not present clearly a break. Only a lower limit of tj > 25 days was proposed by Frail 
et al. (2000). Furthermore, the light curve of its optical afterglow is proportional to i -1 ' 1 , 
in which case no break appears (Galama et al. 1998). We thus include only those bursts 
whose temporal breaks in optical afterglow light curves were well measured in our analysis. 
We obtain a sample of 17 GRBs excluding GRB970508. They are listed in Table 1. 

We correct the observed fluence in a given bandpass to a "bolometric" bandpass of 
1 — 10 4 keV with spectral parameters. The fluence and spectral parameters for a burst fitted 
by different authors may be affected by different criterions (or systematical biases) in their 
works. We thus collect a couple of fluence and spectral parameters from the same original 
literature. For GRB 970828, GRB 980703, GRB 991216, and GRB 011211, the fluence and 
spectral parameters are unavailable in the same original literature, so we choose their fluences 
measured in the widest energy band available in the other literature. 

For GRB 011211, we approximately take the high-energy spectral index f3 to be —2.3 
because it is unreported in Amati (2004). The spectra of HETE- ^-detected GRB020124, 
GRB020813, GRB021004, GRB030226 and XRF030429 are not fitted by the Band function 
but the cutoff power law model. However, it is appropriate that their corresponding "bolo- 
metric" fluences are calculated by the former model with /3 ~ —2.3, avoiding the potential 
systematical bias which is brought by applying different spectral models for one observed 
sample (Barraud et al. 2003). 

The circumburst densities of several bursts in our sample have been obtained from 
broadband modelling of the afterglow emission (e.g., Panaitescu & Kumar 2002). For the 
bursts with unknown n, we assumed n ~ 3 cm -3 as the median value of the distribution of 
the measured densities, together with a constant fractional uncertainty of 80% (Ghirlanda 
et al. 2004a; DLX04). 



3. Cosmological Constraints 

We first fit the E 1 — E p relation for the cosmology of VL M = 0.27, VL\ = 0.73, and 
h = 0.71, and obtain a = 1.53, C = 0.97, a a /a = 0.05 and a c /C = 0.08, together with the 
xl = 21.93/15.0 = 1.46, using Eq. (4) for the estimation of <J E JE 1 (Press et al. 1999). 
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Substituting the best-fit results into Eqs. (7) and (8), we plot the Hubble diagram for the 
observed GRB sample in the concordance model of Qm — 0.27, which is shown in Fig 1 
(filled circles). Comparatively shown are the Hubble diagrams for the binned gold sample 
of SNe la (open circles; Riess et al. 2004), and the theoretical model of VLm = 0.27 and 
Q A = 0.73 (solid line). 

The constraints with Method I are shown in Fig 2 (solid contours). In the concordance 
model of Vt M = 0.27, we find a x 2 do f = 17.91/15 = 1.19. We measure tt M = 0.15±g;^ (la) for 
a flat ACDM universe. DLX04 proposed Method II based on the principle that if there are 
unknown cosmology-independent parameters in the x 2 statistic, they are usually marginal- 
ized over (i.e. integrating the parameters according to their probability distribution). Note 
that DLX04 let the parameter a of the Ghirlanda relation intrinsically equal to 1.50. For 
the purpose of universality, in this paper we let it vary freely, similar to the parameter C of 
this relation (see details in §2.1). The constraints with Method II are shown in Fig 3 (solid 
contours). This method provides a more stringent constraint of 0.14 < Vl M < 0.69 (la) for 
a flat universe. 

By Method I, we recalibrate the Ghirlanda relation for each cosmology, i.e., Q M and 
taken from to 1. We find the half-opening angles of all the bursts are less than 0.23 rad, 
and that the mean a a /a and ac/C are 0.049 and 0.083. As a result, the typical error terms in 
Eq. (8) are asjZS^ ~ 0.052, a k /2k ~ 0.025, WCg/(2 - 2VC5)](3<7 t ./f j ) ~ 0.091, [VU~ e /(2- 
2VC~e)](a n Jn ) ~ 0.118, [1/(2 - 2^C~e)\ (°c/C) ~ 0.055, a[l/(2 - 2^Q,)](a E o bB / E°J> S ) ~ 
0.168, and a a [l/(2 - 2^/Q,)} \ln[E° bs (l + 2)/100]| ~ 0.056. These error terms give a typical 
uncertainty of apparent DM a^ ohs ~ 0.5 magnitudes, which is a factor of ~ 2 larger than 
that derived from the SN la gold sample. 

Firmani et al. (2005) proposed the Bayesian approach to use GRBs as cosmic rulers. 
Being different from their work, we redefined the conditional probability in their method 
and called it Method III in this paper. The constraints on cosmological parameters with this 
method are shown in Fig 4 (solid contours). Compared with previous two methods, Method 
III does give much more stringent confidence intervals at the same confidence levels (C.L.). 
The results are inspiring and demonstrate the advantages of high- ,2 distance indicators in 
constraining cosmological parameters. The data set is consistent with the concordance model 
of tt M = 0.27, yielding a xLf = 19.08/15 = 1.27. We also find 0.16 < tt M < 0.45 (la) for a 
flat universe. 
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4. Simulations and Cosmological Constraints 

4.1. Procedure of Simulations 

As discussed above, the method for the GRB cosmology is different from that for the SN 
cosmology due to the lack of the \ow-z calibration. The \ow-z calibration shall greatly en- 
hance the cosmographic capability of GRBs (assuming no cosmic evolution for the Ghirlanda 
relation). Therefore, one may ask: could a large high- z GRB sample effectively measure cos- 
mology? To answer this question, we carry out a Monte Carlo simulation and use a large 
simulated sample to investigate the cosmographic capabilities in different scenarios. 

Our simulations are based on the Ghirlanda relation which is calibrated for Qm = 0.27 
and £7a = 0.73. Making use of the sample in Table 1, we find a = 1.53, C = 0.97, 6 < 
0.2 rad, and log(£ 7 /l erg) G [49.82,51.96]. We also find \og{E° bs /l KeV) € [1.54,2.89]. 
These restrictive conditions are imposed upon our simulations. Each simulated GRB is 
characterized by a set of S' b , t'j, n', z and E' , where S' b is the bolometric fluence. 

The simulation procedure is as follows: 

(1) We consider the lognormal distributions for the observables Sb, tj and n. From 
the 17 GRBs, we find (log (S^/l erg cm -2 )) ~ —4.46 with cr log (5 6 / lergcm -2) ~ 0.78, and 
(log (tj/1 day)) ~ 0.17 with <7io g (t-/iday) ~ 0.27. Because n is unavailable in the literature for 
most GRBs, we take (log (n/1 cm -3 )) ~ 0.40 with o"i og (n/i C m- 3 ) ~ 0.25. The observable z is 
selected according to its observational distribution cut by an upper limit z ~ 4.5, which is 
shown in Fig 5. The uncertainty of z is ignored. 

(2) We also consider the lognormal distributions for the fractional uncertainties of the 
observables Sb, tj, n and E p . From the observed sample, we find (log (o~s b /Sb)) ~ —0.97 
with ai og[aSb/Sb } ~ 0.16, (log (<J t Jtj)) ~ -0.91 with a log[at ^ /tj] ~ 0.45, (log (<7 n /n)) ~ -0.30 
with 0~\ O g[ an / n ] ~ 0.10, and (log (a Ep /E p )) ~ -0.83 with 

^QgWEp/Ep] ~ 0.26, respectively. We 
ensure in code that the fractional uncertainties of the observables Sb, tj, n and E p are less 
than 25%, 35%, 100% and 35%, respectively. 

(3) We simulate a GRB characterized by a set of (Sb ± cg b , tj ± a tj , n ± o n and z) 
according to the distributions that these parameters follow, compute its E 1 in the cosmic 
model of Qm = 0.27 and = 0.73, and then calculate its E p by the Ghirlanda relation of 
(£ 7 /10 50 ergs) = 0.97(£ p /100keV) L53 . 

(4) The GRB generated in step (3) follows the "rigid" Ghirlanda relation. We add a 
random deviation to each parameter, except for z, to make this burst more realistic, i.e., S' b = 
S b + l.l(-l) m a Sb , t'j = tj + l.l(-l) m a tj , n' = n + l.l(-l)"V n and E' p = E p + l.l(-l) m a Ep , 
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where m is randomly taken from and 1. 

(5) Using the parameters S' b , t'j, n', z, and E' p , we compute the parameters 9' and E' 
for tt M = 0.27 and fi A = 0.73, and the quantity E' p obs = E' p /(1 + z). 

(6) Since 9 < 0.2 rad, log(£ 7 /l erg) G [49.82,51.96] and log(E° bs /l KeV) G [1.54,2.89] 
are valid for the observed sample, we require that 9', log (E' /l erg), and \og(E' p obs /l KeV) of 
a simulated GRB must be within the corresponding ranges. 

(7) Repeat steps (3)- (6) to generate a sample of 80 bursts. 

The half-opening angles of the simulated sample are less than 0.23 rad when Qm and 
Q\ are taken from to 1, as in the observed sample. We carry out a circular operation to 
achieve the simulated sample, which may be established in the Swift era. 

4.2. Constraints from Simulated GRBs 

By Method I, we find the typical a a /a and a c /C for the large sample decrease to 0.02 
and 0.05. This is because for the small observed sample, the dispersion of the Ghirlanda 
relation is mainly contributed by a few "outliers", while for the large simulated sample, the 
bursts are distributed around the "rigid" Ghirlanda relation with a Gaussian distribution 
(see §4.1). Such a large sample seems to be more realistic. 

Among several scenarios, we only perform to what degree the simulated sample can 
constrain the Qm — ^a parameters with Method III. The results are shown by solid confidence 
intervals in Fig 6. Also shown are the constraints derived from the SN gold sample (dashed 
contours) at the same confidence levels. The main points revealed by this figure are: (1) 
A large high-z GRB sample could effectively measure cosmology. The simulated GRB data 
are consistent with the concordance model of VL M ps 0.3, yielding a xLf — 94. /78. 1.20. 
With the prior of a flat universe, the mass density is tt M = 0.30i{J;{jjj at the 68.3% level. (2) 
GRBs can well constrain Q, M rather than due to their high redshifts. The orientation 
of the elliptical contours is almost vertical to the Qm axis. This is an advantage of GRBs 
over SNe la for the cosmological use. (3) GRBs supply complemental content to the SN 
cosmology. The 157 SNe la provide evidence of cosmic acceleration at a very high confidence 
level. A large GRB sample would reach a similar conclusion. Alone the SN sample nearly 
rules out the flat universe model at la level, but a combination of GRBs and SNe makes the 
concordance model of Qm ~ 0.3 more favored and thus more agreement with the conclusion 
from WMAP observations (Spergel et al. 2003). 

If the \ow-z calibration is realized for GRBs, then the solid confidence regions in Fig 
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6 will become smaller so that they lie in the part of cosmic acceleration at a high C.L. in 
the Qm — plane (not shown in this work). We here present a rough estimation of the 
detection rate of \ow-z bursts in the Swift era. Taking Qm — 0.27, fl\ — 0.73, and h = 0.71, 
the observed rate of bursts with redshift less than z is 



dN r 

~dt~ Jo 



dz dV(z)R GRB (z) 
dz 1 + z 



where Rgrb(z) is the comoving GRB rate density, and dV(z)/dz is the isotropic comoving 
volume element, 



dV A 

— = 4:71 

dz 



dz' 



y/n M (l + Z'f + Qf 



(18) 



y/n M (i + z) 3 + n A ' 



We assume Rgrb(z) oc Rsn(z) oc Rsf(z), where Rsn(z) and Rsf(z) are the comoving rate 
densities of core-collapse supernovae and star formation, respectively (Porciani & Madau 
2001). We also assume the constant fraction k = RqrbI Rsn ~ 10~ 5 , which considers 
the GRB formation efficiency out of the core-collapse SNe and the GRB beaming effect. 
Additionally, we take Rsn(z) — 0.01Rsf(z)M.q. The global star formation rate for the 
Einstein-de Sitter universe (Steidel et al. 1999) is 

Thus, we derive the detective probability of low- z GRBs in the Swift era (2—4 years): when 
z < 0.1, there will be < 1 burst; and when z < 0.2, there will be a few bursts. These results 
agree with the present observational GRB data. In addition, taking into account possible 
different origins of very low- z bursts (e.g. GRB980425) and high- ,2 bursts, it might not be 
valid to directly apply the low- z calibrated relation to a high- z sample. However, if the 
cosmic evolution of the "standard candle" relation is ignored, then in the future one could 
consider a sample of GRBs with redshift z < 0.2 for the \ow-z calibration. As the zero-order 
approximation, the theoretical luminosity distance of such a GRB sample can be written as 
d L (z) = z(c/H ). 



5. Conclusions and Discussion 

At present, GRBs with known redshifts are about 45 out of a few thousands, among 
which 17 are available to derive the Ghirlanda relation. This relation is so tight that it has 
been considered for the cosmological use. Although the \ow-z calibration of the Ghirlanda 
relation is not realized, the observed 17 high- ,2 bursts still provide interesting and even 
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inspiring results. We find GRBs independently place a constraint on the mass density 0.16 < 
tt M < 0.45 (la) for a flat ACDM model. 

GRBs are becoming more and more standardized candles. Through the Ghirlanda 
relation, the mean scatter in GRBs is a factor of ~ 2 larger than that in SNe la. However, the 
disadvantage of larger scatter in GRBs is compensated, to some extent, by their advantages 
of high redshifts and immunity to the dust extinction. The shape of the constraints in 
Fig 6 implies that GRBs could not only measure the mass density Qm but also provide 
complementarity to the SN cosmology. The Swift satellite is hopeful to establish a large 
GRB sample with known redshifts, perhaps including \ow-z bursts. In this sense, the GRB 
cosmology now lies in its babyhood. 

A reliable theoretical basis of the Ghirlanda relation is also important for GRBs as a 
cosmic ruler. Using the small angle approximation for all bursts, a scaling analysis gives 
E 1 oc Ep (a ~ 1.5) as long as the parameters such as the spectral index p of the distribution 
of accelerated electrons, the the energy equipartition factor e e of the electrons, the energy 
equipartition factor eb of magnetic field, the bulk Lorentz factor 7, etc., or their combinations 
are clustered (Zhang & Meszaros 2002; Dai & Lu 2002; Wu, Dai & Liang 2004). This 
power-law relation could also be understood by the emissions from off-axis relativistic jets 
(Yamazaki, Ioka & Nakamura 2004; Eichler & Levinson 2004; Levinson & Eichler 2005) and 
the dissipative photosphere model (Rees & Meszaros 2005). Different plausible explanations 
imply that this topic needs further investigations. 

It should be pointed out that the Ghirlanda relation is obtained under the framework 
of the uniform top-hat jet model. Other input assumptions include the uniform circumburst 
medium density n and the constant efficiency r^ 7 of converting the initial ejecta's kinetic 
energy into 7-ray energy release. However, n and i] y should be different from burst to burst, 
and n is variable for a burst in the wind environment (Dai & Lu 1998; Chevalier & Li 1999). 
Thus, the quantity n?7 7 in the Ghirlanda relation might not be clustered for the observed 
bursts (Friedman & Bloom 2004). New relations with as few well-observed quantities as 
possible are required for improvement 2 . In this paper, for those bursts with unknown n, 
we assumed n ~ 3 cm -3 as the median value of the distribution of the measured densities, 
together with a constant fractional uncertainty of 80% (Ghirlanda et al. 2004a; DLX04). 
We also treat r]~ ( = 0.2 and a Vl = for all the 17 bursts. 

Finally, the Ghirlanda relation is not valid for all the observed BATSE sample. This 



2 After the submission of this paper, Liang & Zhang (2005) and Xu (2005) proposed new relations between 
the isotropic 7-ray energy and the vF v peak energy by considering the break time of an afterglow light curve. 
Clearly, these three quantities are directly observed. 



-14- 



implies that this relation may suffer from the data selection effect (Band & Preece 2005). Its 
validity will be tested by the ongoing observations of the Swift mission. However, no matter 
whether the low-z sample is established or not, the GRB cosmology is expected to progress 
much in the coming years. 
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a The spectral parameters fitted by the Band function. The fractional uncertainties of E° bs are taken as 20% when not reported, and the fractional uncertainty 
of fc-correction is fixed as 5%. 

b The fluences and their errors in the observed energy band. The fractional errors are taken as 10% when not reported. The fluence and spectral parameters 
of a GRB are selected from the same original literature as possible. If this criterion unsatisfied, fluence are chosen in the widest energy band. 

c Aftcrglow break times and errors in the optical band. 

d The circumburst densities and errors from broadband modelling of the afterglow light curves. If no available the value of n is taken as 3.0 ± 2.4 cm 
e The constant efficiency »7 7 of converting explosion energy into 7-ray emission for each burst. 
f References in order for 2, E° ba , [a,/3], S~ f , tj, and n. 

References. — (1) Djorgovski et al. 2001; (2) Djorgovski ct al. 1998; (3) Kulkarni ct al. 1999; (4) Vreeswijk et al. 2001; (5) Lc Floc'h ct al. 2002; (6) 
Vreeswijk et al. 2001; (7) Piro ct al. 2000; (8) Holland ct al. 2002; (9) Hjorth et al. 2003; (10) Price et al. 2003a; (11) Barth ct al. 2003; (12) Matheson et al. 
2003; (13) Vreeswijk et al. 2003; (14) Grciner ct al. 2003a; (15) Rol ct al. 2003; (16) Greiner et al. 2003b; (17) Jakobsson, et al. 2004; (18) Jimenez et al. 2001; 



(19) Amati et al. 2002; (20) Amati 2004; (21) Barraud et al. 2003; (22) Price et al. 2003a; (23) Sakamoto et al. 2004b; (24) Crew et al. 2003; (25) Vanderspek 
ct al. 2004; (26) Bloom et al .2003; (27) Frail et al. 2003; (28) Kulkarni et al. 1999; (29) Stanek ct al. 1999; (30) Masctti et al. 2000; (31) Bjornsson ct al. 2001; 
(32) Halpern ct al. 2000; (33) Jakobsson ct al. 2003; (34) Berger et al. 2002; (35) Barth ct al. 2003; (36) Holland et al. 2003; (37) Holland ct al. 2004; (38) 
Klosc ct al. 2004; (39) Andersen ct al. 2003; (40) Price ct al. 2003b; (41) Panaitcscu & Kumar 2002; (42) Schaefer et al. 2003; (43) Ticngo ct al. 2003 
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Fig. 1. — Hubble diagrams for the observed GRB sample (filled circles) and for the binned 
gold sample of SNe la (open circles). The GRB E 7 — E p relation has been calibrated in the 
cosmic model of Qm = 0.27, fl^ = 0.73 and h = 0.71 (solid line). 
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Fig. 2. — Joint confidence intervals (68.3%, 90% and 99%) in the Vt M — f2 A plane from the 
17 GRBs with Method I. 
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Fig. 3. — Joint confidence intervals (68.3%, 90% and 99%) in the Vt M — plane from the 
17 GRBs with Method II. 
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Fig. 4. — Joint confidence intervals (68.3%, 90% and 99%) in the Vt M — plane from the 
17 GRBs with Method 111. 
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Fig. 5. — Histogram for the redshifts of 80 simulated GRBs, following the redshift distribu- 
tion of the GRBs observed sofar. 
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Fig. 6. — Joint confidence intervals (68.3%, 95.4% and 99.7%) in the Ojh-Oa plane from the 
80 simulated GRBs in this work (solid contours) and from the 157 SNe la in Riess et al. 
(2004) (dashed contours). The black dot marks VL M = 0.30 and £1\ = 0.70. 



